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ABSTRACT 

Context. Optical interferometry allows a measurement of the intensity profile across a stellar disc, leading to a direct test and calibration of 
theoretical model atmospheres as well as to a precise determination of fundamental stellar parameters. 

Aims. We present a comparison of the visual and near-infrared intensity profile of the M0 giant y Sagittae to plane-parallel ATLAS 9 as well as 
to plane-parallel & spherical PHOENIX model atmospheres. 

Methods. We use previously described visual interferometric data obtained with the Navy Prototype Optical Interferometer (NPOI) in July 
2000. We apply the recently developed technique of coherent integration, and thereby obtain visibility data of more spectral channels (526- 
852 nm) and with higher precision than before. In addition, we employ new measurements of the near-infrared K-band (~ 2200 nm) diameter 
of y Sagittae obtained with the instrument VINCI at the ESO VLT Interferometer (VLTI) in 2002. 

Results. The spherical PHOENIX model leads to a precise definition of the Rosseland angular diameter and a consistent high-precision diameter 
value for our NPOI and VLTI/VINCI data sets of 0r oss = 6.06 ± 0.02 mas, with the Hipparcos parallax corresponding to Rr oss = 55 ± 4 R B , 
and with the bolometric flux corresponding to an effective temperature = 3805 ± 55 K. Our visual visibility data close to the first minimum 
and in the second lobe constrain the limb-darkening effect and are generally consistent with the model atmosphere predictions. The visual 
closure phases exhibit a smooth transition between and n. 

Conclusions. The agreement between the NPOI and VINCI diameter values increases the confidence in the model atmosphere predictions from 
optical to near-infrared wavelengths as well as in the calibration and accuracy of both interferometric facilities. The consistent night-by-night 
diameter values of VINCI give additional confidence in the given uncertainties. The closure phases suggest a slight deviation from circular 
symmetry, which may be due to surface features, an asymmetric extended layer, or a faint unknown companion. 

Key words. Techniques: interferometric - Stars: late-type - Stars: AGB and post-AGB - Stars: atmospheres - Stars: fundamental parameters 
- Stars: individual: y Sagittae 



1. Introduction 

Cool giants on the red giant branch (RGB) and asymptotic 
giant branch (AGB) are very luminous and extended, have a 
low surface temperature, and their atmospheres can thus be 
rich in molecules. Cool giants are the most important source 
of dust formation and its delivery to the interstellar medium. 
The detailed structure of their extended atmospheres, includ- 

* Based on data obtained with the Navy Prototype Optical 
Interferometer (NPOI). The NPOI is a joint project of the Naval 
Research Laboratory and the United States Naval Observatory in co- 
operation with Lowell Observatory, and is funded by the Office of 
Naval Research and the Oceanographer of the Navy. 
** Based on public commissioning data released by the Paranal 
Observatory, Chile, and received via the ESO/ST-ECF Science 
Archive Facility. 



ing the effects from circumstellar molecular and dust layers, 
are still a matter of investigation and debate (cf., e.g. Scholz 
IT985l[T^8ll2'no'T1 Perrin et al. EDTHI Ohnaka llUOlal Ireland & 
Scholz l2006l . 

Theoretical atmosphere models predict in general the spec- 
trum emerging from every point of a stellar disc. Optical in- 
terferometry provides the strongest observational constraint of 
this prediction by resolving the stellar disc. In addition, the con- 
straints on the intensity profiles allow us to find meaningful 
definitions of the stellar radius and its precise measurement. 

For regular cool non-pulsating giants, the centre-to-limb 
variation (CLV) is mainly characterised by the limb-darkening 
effect, which is an effect of the vertical temperature profile 
of the stellar atmosphere. The strength of the limb-darkening 
can be probed by optical interferometry in two ways (cf., e.g. 
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Hanbury Brown et al. 119741 Quirrenbach et al. 119961 Burns 
et al. nW71 Hajian et al. [TWHl Wittkowski et al. I200T1 120051 
Aufdenberg et al. 2005): (1) by measuring variations of an 
equivalent uniform disc diameter (i.e. the uniform disc that has 
the same integral flux as the true intensity profile) as a func- 
tion of wavelength, and (2) by directly constraining the star's 
intensity profile in the second and higher lobes of the visibility 
function at one or several bandpasses. 

It was found that pulsating giants as well as supergiants 
may exhibit more complex intensity profiles at near- and mid- 
infrared wavelengths, showing Gaussian-shaped intensity pro- 
files, tail-like extensions to a photospheric intensity profile, and 
multiple components, such as a photosphere plus a circumstel- 
lar shell (cf., e.g. Woodruff et al. 120051 Ohnaka 2004al l2004bl 
Perrin et al. 120051 120051 Fedele et al. 120051 . Additionally, ob- 
served intensity profiles might be affected by dust shells (e.g. 
Ohnaka et al. 2005 Ireland & Scholz 2006 1 or horizontal sur- 
face inhomogeneities (e.g. Burns et al. 1997). 

In Wittkowski et al. (2001 hereafter Paper I), we used the 
Navy Prototype Optical Interferometer (NPOI, Armstrong et al. 
I1998> . used the method of baseline bootstrapping (cf. Hajian et 
al. Q998 1, and developed improved methods of compensation 
of noise and detection bias terms, in order to obtain precise vi- 
sual visibility measurements in the second lobe of the visibility 
function for three cool giants. We found agreement with pre- 
dictions by plane-parallel ATLAS 9 ( 1993 ) model atmospheres 
within the obtained wavelength range and precision. Thereby, 
the strength of the limb-darkening effect and the stars' funda- 
mental parameters were constrained. Aufdenberg & Hauschildt 
(2003 1 compared one of the NPOI observations of y Sagittae 
from Paper I to a spherical PHOENIX ( 1999 1 model atmosphere 
and found agreement. In Wittkowski et al. ( 120041 hereafter 
Paper II), we directly measured the limb-darkening effect of 
the M4 giant \j/ Phoenicis using the ESO Very Large Telescope 
Interferometer (VLTI) in the near-infrared /f-band, confronted 
the observations with predictions by independently constructed 
ATLAS 9 and PHOENIX model atmospheres, and found agree- 
ment with all considered models. 

Recently, Hummel et al. (2003) developed the method of 
coherent integration and its application to NPOI data in order to 
increase the precision of visibility measurements. This method 
was recently applied by Peterson et al. ( 2006a 2006b I to NPOI 
observations of Altair and Vega. 

Here, we reanalyse the NPOI data of the M0 giant 
y Sagittae (HR7635, HD 189319), the brightest of the targets 
in Paper I, using the newly developed method of coherent inte- 
gration. We obtain visibility data with higher precision than in 
Paper I, and -due to the lower noise- are able to make use of 
more spectral channels toward the blue end of NPOI's wave- 
length range. Now, the wider wavelength range covers 526- 
852 nm, compared to 649-852 nm in Paper I. We thus also 
increase our maximum spatial resolution from «3.3mas to 
«2.7mas, which gives important additional visibility data in 
the second lobe that are sensitive to the limb-darkening ef- 
fect. In addition, we observed y Sagittae with the ESO Very 
Large Telescope Interferometer (VLTI) and its /^-band instru- 
ment VINCI, in order to compare results derived from differ- 
ent interferometric facilities, and to probe the consistency of 



the wavelength-independent Rosseland diameter from visual to 
near-infrared wavelengths. 

The cool giant y Sagittae does not appear in the Combined 
General Catalogue of Variable Stars (Samus et al. 2004), in- 
dicating that it lacks strong photometric variability. Thus, it 
is a good target for the purpose of calibrating model atmo- 
spheres and deriving high-precision fundamental parameters. 
The spectral type has been listed as K5-M0III by Morgan 
& Keenan ( 119731 ). and been revised to M0 III by Keenan & 
McNeil ( 1989). Wisniewski & Morrison (private communica- 
tion) confirm by means of optical echelle spectra recently ob- 
tained at Ritter Observatory that ySagittae's spectrum closely 
resembles that of the MK standard /iUMa (M0 III). We de- 
termine the bolometric flux of y Sagittae to /t, i = (2.57 + 
0.13) x 10~ 9 W/m 2 by means of a spline fit and integration of 
the narrow-band spectrophotometric data by Alekseeva et al. 
( 1997) covering 405 nm to 1080 nm complemented by broad- 
band photometry shortward and longward of this range from 
the 13 -colour photometry by Johnson et al. dl975> . The val- 
ues for / bol of (2.79 + 0.14) x 10"^ W/m 2 and (2.83 + 0.14) x 
10~ 9 W/m 2 by Alonso et al. (1999) and Mozurkewich et al. 
(2003), respectively, are derived from broad-band photome- 
try alone and likely overestimate /t, i because of a too sparse 
sampling of the visual spectrum including the TiO band heads 
and other features. The limb-darkened angular diameter of 
y Sagittae has been determined in Paper I to be 6. 18 +0.07 mas, 
based on a comparison of NPOI visibility data to ATLAS 9 
model atmospheres. This value corresponds to a limb-darkened 
radius of 56+4 R Q , derived with the Hipparcos parallax of 
1 1 .90 ± 0.7 1 mas (Perryman & ESA. 1T9971 . These values of an- 
gular diameter, absolute radius, and bolometric flux constrain 
the effective temperature to T e s = 3768 K ± 70 K, and the lu- 
minosity to log L/L e = 2.75 ± 0.10. Placing y Sagittae on the 
Hertzsprung Russel diagram using these values, and compar- 
ing to stellar evolutionary tracks by Girardi et al. (|2000 ) as in 
Paper II (Fig. 1 of Paper II) we can estimate a mass of M — 
1.3 + 0.4 M , and thus a surface gravity of log g=1.06 + 0.22. 
These values are used as an a-priori estimate for our analysis 
and will be refined in the conclusions. 

2. NPOI measurements 

2.1. NPOI observations 

We reanalyse the visual y Sagittae data in PapeiT obtained with 
NPOI on July 21, 2000. The centre (C), east (E) and west (W) 
siderostats of the astrometric sub-array of NPOI were used to 
obtain baselines of ground length 18.9 m (CE), 22.2 m (CW), 
and 37.5 m (EW). The data were recorded in 32 spectral chan- 
nels of equal width in wavenumber and covering the band from 
« 450 nm to 850 nm. Due to low photon count rates only the 
10 reddest channels could be used in Paper I (covering 649 nm 
to 852 nm). We reanalyse the same raw data using the newly 
developed coherent integration algorithm as first described by 
Hummel et al. (2003). The details of our new analysis are de- 
scribed below. The benefits of the new analysis include an im- 
proved signal-to-noise-ratio (SNR) of the visibility data on the 
long EW baseline, as well as a much improved SNR of the 
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Fig. 1. Incoherent integration. Off-fringe squared visibility am- 
plitude, i.e. the visibility bias that remains for data off the fringe 
packet and that is compensated after the Z 2 compensation (see 
text for more details). This bias can be described by a power 
law as a function of photon rate P. As an example, the residual 
bias is shown for the four reddest channels on the EW baseline. 
Data are 2 ms incoherent integrations from July 22. Channels 1 
though 4 use symbols plus, star, diamond, and triangle, respec- 
tively. Power-law fit coefficients are given for each channel, in 
the same order as shown on the corresponding plot for the co- 
herent analysis in Fig.|2] 

triple amplitudes and phases. These improvements enable us to 
use the 20 reddest channels in the present paper, now covering 
526 nm to 852 nm. 

2.2. NPOI data reduction and calibration 

Coherent integration The NPOI detector configuration used 
for our data set (July 2000) sampled a single fringe of each 
two telescope interference pattern using 8 bins every 2 ms. 
This time interval is called the instrumental coherent (because 
it is phase preserving) integration time. Increasing this time 
would eventually lead to a complete loss of fringe contrast due 
to atmospheric fringe motion which is not perfectly compen- 
sated by an interferometer such as NPOI using the group delay 
method for fringe tracking. 

Therefore, in the so-called incoherent analysis, as used in 
Paper I, the bin counts of each 2 ms sample would be Fourier 
transformed, and an unbiased estimate for the squared visibility 
amplitude derived as follows (see Shao et al. 1988 1: 



W\ 2 



< X 2 + Y 2 - o-){N) > 
<N> 2 ' 



(1) 



where X and Y are the real and imaginary parts of the Fourier 
transform, respectively, and cr 2 is the variance of the inten- 
sity caused by photon and detection noise. In the case of pure 
Poisson noise cr 2 equals N, the total number of bin counts of 
each sample. In the case of NPOI, the data of which show non- 
Poisson noise due to afterpulsing of the APDs (avalanche photo 



Fig. 2. Coherent integration. As Fig. [2 but for 200ms coherent 
integrations as used in this paper. For reasons of comparison, 
the bin counts are renormalised to 2 ms intervals. It is an addi- 
tional benefit of the coherent integration that this residual bias 
shown here is clearly reduced compared to the incoherent inte- 
gration in Fig.n 

diodes), the bias is estimated according to procedures described 
in PapeiT (see also below the paragraph "Correction of noise 
and detection bias terms"). 

The signal to noise ratio SNR of the squared visibility esti- 
mator is as follows (Shao et al. 1988 1: 



SNR(V 2 ) = ^M l/2 NV 2 



1 ? 
1 + -NV 2 
2 



-1/2 



(2) 



where M is the number of samples averaged. One can see that 
by increasing the coherent integration time, i.e. increasing N, 
rather than increasing M a larger gain in SNR can be realised. 
This is true as long as NV 2 is much smaller than unity, oth- 
erwise nothing can be gained by a coherent average over an 
incoherent average. Here, we chose a coherent integration time 
of 200 ms which still results in NV 2 > 5 for photon rates of 10 
per 2 ms at squared visibilities of about 0.005. Subsequently, 
the squared visibility was computed for the coherent samples, 
and then incoherently averaged in 2 s intervals. The complex 
triple products were computed from the coherent samples as 
well, but vector averaged to preserve the phase. Coherent in- 
tegration for time intervals longer than the instrumental coher- 
ent integration time (we chose 200 ms) require the alignment 
of the complex visibilities, equivalent to removing the relative 
offsets between successive samples of the fringe. As described 
by Hummel et al. (2003 (, this can be done in the off-line data 
analysis by making use of the visibility phase derivative with 
wavenumber. This quantity, the so-called group delay, is zero 
for the white light fringe which is located at zero relative opti- 
cal path length difference where the fringes of all colours inter- 
fere constructively. 

For the observations described here on y Sagittae, the im- 
portance of coherent integration follows from the very small 
visibility amplitudes measured on the long EW baseline since 
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it samples the second lobe of the Fourier transform of the stellar 
disc brightness profile. This baseline is therefore most sensitive 
to stellar limb darkening, the focus of this work. While the low 
visibility amplitudes on this baseline would prevent a precise 
determination of the group delay, the following paragraph de- 
scribes how to obtain this estimate in a different way. 

Phase bootstrapping A design feature of the NPOI array 
(Armstrong et al. 1998) is the ability to realise configurations 
which allow long baselines to be boot-strapped by shorter base- 
lines. By this we mean the ability to track and observe fringes 
on long baselines even though the fringe contrast can be too low 
to allow detection of the fringes for tracking purposes. This is 
achieved by detecting and tracking fringes on the shorter base- 
lines, and making use of the fact that the sum of all fringe de- 
lays over baselines in a closed loop must equal zero. 

In the case described here, the EW baseline is boot-strapped 
by the CE and CW baselines. Therefore, the fringe delay on 
the EW baseline is equal to the difference of the delay between 
the other two baselines, and can thus be computed in this way 
without using the measurement on the EW baseline itself. 

Correction of noise and detection bias terms As shown in 
Paper I, the standard correction for non-Poisson statistics of the 
NPOI detectors used by Hummel et al. i 1998b (for a] of Eq.[T] 
useZ 2 withZ = B x -B 2 +Bj,-B A +B 5 -B( s +B 1 -B ? „ where 5, are 
the bin counts) can be improved by compensating a small resid- 
ual bias that remains after the Z 2 compensation. This residual 
bias is calibrated by a power law as a function of counts by ob- 
serving off the fringe for a number of stars of different bright- 
ness (cf. Eqs. 3 and 4 of Paper I). This additional bias correction 
is important because of the small visibility amplitudes mea- 
sured on resolved stars, where a constant residual bias would 
strongly contribute to the visibility values. A welcome side ef- 
fect of the coherent integration is that, because the number of 
photons counted in a coherent integration, N, increases, the ra- 
tio of the bias to the signal decreases. Therefore, at a count rate 
high enough, the remaining bias after the Z 2 compensation be- 
comes negligible. Since we chose a coherent integration time 
of 200 ms, a measurable, though much smaller, residual bias 
than in PapeiT remained. This residual bias was removed using 
exactly the same procedures as described in Sect. 3 of PapeiT 
by reducing data obtained on July 22 during which off-fringe 
data were recorded. Figures ^ an d Q show as an example for 
the 4 reddest channels of the EW baseline the remaining bias 
after the Z 2 compensation, obtained for the incoherent analysis 
as used in Paper I (2 ms integrations) and for the coherent anal- 
ysis used here (200 ms coherent integrations, and normalised 
to a 2 ms interval), respectively. It can be seen that the resid- 
ual bias for the coherent analysis is clearly reduced compared 
to the incoherent average. This is an additional benefit of the 
method of coherent integration. 

Calibration The calibration of the data followed the same pro- 
cedures as used in Paper I. The B6 giant eDelphini, 12 de- 
grees away and observed in an interleaved way during the same 



night, is used as the main calibration star. The diameter of 
eDelphini is estimated to be 0.3 mas based on a calibration of 
the visual magnitude and (R - 7) colour index by Mozurkewich 
et al. (I1991> . As secondary calibrators we used i 2 Cygni and 
7r 2 Pegasi, both 32 degrees away from y Sagittae. The calibra- 
tion errors of the squared visibility amplitudes (derived from 
the scatter of the primary and secondary calibrator amplitudes 
and added in quadrature to the formal, i.e. photon noise induced 
errors) are now larger than those quoted in Paper I for the boot- 
strapped EW baseline, namely 13% versus 7%, while the cal- 
ibration error for the other two baselines was only marginally 
larger. This is acceptable since the statistical errors dominate 
the total error budget for the EW baseline. 

2.3. NPOI results 

As a first characterisation of our NPOI data, we use models of 
a uniform disc (UD, 7=1 for < // < 1,7 = otherwise), and 
a fully darkened disc (FDD, 7 = jj). Here, 7 is the intensity, fi = 
cos (or n = y/l - JrjRp) the cosine of the angle between the 
line of sight and the normal of the surface element of the star 
(R the stellar radius, r the distance from the centre of the disc). 
Monochromatic synthetic visibility values V were obtained for 
the UD and FDD cases, and subsequently integrated over the 
bandpass of each NPOI spectral channel (covering frequencies 
V] to v%) as 



V t = 



X; F v V(v)dv 



(3) 



with F Y the flux from the blackbody radiation. We use = 
3768 K (see Sect. [I). Variations of T e s within its uncertain- 
ties do not have a significant effect on the visibility values. 
Finally, the synthetic squared visibility values for each of the 
baselines (\Vi(CE)\ 2 , |V,(CW)| 2 , |V ; (7iW)| 2 ), the triple ampli- 
tude (Vf EW = \Vi(CE) Vi(CW) Vi(EW)\), and the closure phase 
(Of EW = if Vi(CE) VdCW) Vi(EW) > 0, <S>f EW = n if 
Vi(CE) Vi(CW) Vi(EW) < 0) were obtained. 

Note that in the case of NPOI the monochromatic visibility 
amplitudes are integrated before building the square, while for 
VLTI/VINCI the monochromatic squared visibility amplitudes 
are integrated (cf. PapeiTI). The reason for the difference is 
that the data processing of NPOI first integrates the photons 
on the APDs and the squared visibility is computed from the 
already integrated bin counts, while for VLTI/VINCI first the 
full powerspectrum is computed and integrated thereafter. This 
leads to noticeable differences, in particular around the minima 
of the visibility function. 

Best fitting angular diameters Oud.fdd are derived from a 
least square optimisation. The resulting values are shown in 
Table n together with the reduced^ 2 values. The number of 
degrees of freedom is 644 (7 observations times 19 spectral 
channels times (3 squared visibility plus 1 triple amplitude plus 
1 closure phase), minus 21 values flagged for quality). The 
formal errors of the obtained diameter values are of the order 
of 0.01 mas, and are small compared to calibration uncertain- 
ties that are estimated to ~ l%~0.06mas. Total errors are thus 
~ 0.06 mas. 
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Fig. 3. Squared visibility amplitudes of ySge obtained from NPOI on the EW baseline. Also shown are synthetic visibility 
curves of the best fitting model atmospheres as described below in Sect.@] Each model fit is performed to all NPOI visibility data 
(squared visibility amplitudes, triple amplitudes, and closure phases) simultaneously. The parameters of the plotted model curves 
are listed in Table[5]. The synthetic visibility values are calculated for each specific bandpass of our observation and these points 
are connected by straight lines. A 7th used NPOI observation is not included for the sake of clarity of the figure and because its 
projected baseline length is much shorter and this observation thus contains little information in the 2nd lobe. 
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Fig. 4. As Fig. [5] but showing the squared visibility amplitudes on the NPOI CW baseline. 



Table 1. Fit results of our NPOI data to models of a uniform 
disc (UD) and a fully darkened (FDD) disc. The formal errors 
of the diameter values are ~ 0.01 mas, additional calibration 
uncertainties are ~ 0.06 mas, total errors thus ~ 0.06 mas. 



Model 




Diameter 


Parameter 


a 




UD 


©UD 


= 5.64 mas 


a = 





11.0 


FDD 


©FDD 


= 6.59 mas 


a = 


1 


5.6 



Figures [3] EJH] [6] and show the obtained NPOI squared 
visibility amplitudes on baselines EW, CE, CW, the NPOI triple 
amplitudes, and the NPOI closure phases, respectively. Also 
shown are the model atmosphere predictions as described be- 
low in Sect.|U 

The gain with respect to Paper I in the signal-to-noise ratio 
and in the number of usable spectral channels is thanks to the 
method of coherent integration, as can be seen by comparing 
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Fig. 5. As Fig. [5] but showing the squared visibility amplitudes on the NPOI CE baseline. 
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Fig. 6. As Fig. 13 but showing the NPOI triple amplitudes. 



these Figs, to the results based on incoherent averaging (Fig. 3 3. VLTI/VINCI measurements 

of Paper I). 

3.1. VLTI/VINCI observations 



The results (Figs. I3f7l and Tabled show that the visual in- 
tensity profile of ySge is limb-darkened, clearly closer to a 
FDD model than to a UD model, while both of these simple 
descriptions do not provide a very good representation of our 
data. A detailed comparison of our visibility data to model at- 
mosphere predictions is discussed below in Sect.|4] 



The near-infrared K-band interferometric data of y Sagittae 
were obtained with the ESO Very Large Telescope 
Interferometer (VLTI, Glindemann et al. 2003 1, the instrument 
VINCI (Kervella et al. 2003 1, and the two VLTI test siderostats 
on June 28, July 8, July 11, July 15, August 8, September 12, 
and September 18, 2002. These data are public commission- 
ing data released from the VLTI 1 . The VLTI stations E0 and 
Gl forming a ground baseline length of 66 m were used for 



http://www.eso.org/projects/vlti/instru/vinci/vinci_data_sets.html 
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Fig. 7. As Fig. |3 but showing the NPOI closure phases. Note that the slope of the model flip from to n is an artifact because 
the model values are only calculated for each NPOI spectral channel. 



Table 2. Characteristics of the stars that were used as cali- 
bration stars for our VLTI/VINCI observations of y Sagittae. 
Listed are the spectral type, the A'-band magnitude, the 
uniform-disc diameter and its error, and the effective temper- 
ature, all from Borde et al. (I2UU2I based on Cohen et al. 119991 . 



Star 


Sp. Type 


K 


©UD 


o-(0) 




56 Aql 


K5 III 


1.76 


2.45 


0.028 


4046 


31 Ori 


K5 III 


0.90 


3.56 


0.057 


4046 


58Hya 


K2.5 Illb 


1.13 


3.13 


0.035 


4318 


66 Aql 


K5 III 


1.76 


2.37 


0.030 


4046 


70 Aql 


K5 II 


1.21 


3.18 


0.037 


4064 


eCen 


K0- Illb 


-0.26 


5.32 


0.058 


4656 


<p l Aqr 


Kl- III 


1.79 


2.18 


0.025 


4508 


AGm 


K3III 


1.44 


2.64 


0.030 


4256 


^Sgr 


Kl Illb 


0.40 


4.13 


0.047 


4508 


n 2 Ori 


K0 Illb 


1.69 


2.14 


0.023 


4656 


^Phe 


K5 III 


1.52 


2.69 


0.032 


4046 



all our observations. The observations were repeated during 
7 different nights spread over more than 2 months in order to 
compute the night-to-night variation of the obtained diameter 
and thereby to estimate the calibration uncertainty caused by 
different atmospheric and possibly instrumental conditions. All 
data were obtained as series of typically 100 or 500 interfero- 
grams with a scan length of 250 fim and a fringe frequency of 
295 Hz. 

The stars 56 Aquilae and 3 1 Orionis were used as primary 
calibration stars and were observed in each of our observation 
nights close in time to the y Sagittae observations. A number of 
additional calibration stars observed during these nights were 
used as secondary calibrators for y Sagittae. The characteristics 
of all calibration stars used are taken from Borde et al. ( 2002 
based on Cohen et al. 1999) and are listed in Table|2] 



3.2. VLTI/VINCI data reduction and calibration 

We computed mean coherence factors for each series of inter- 
ferograms using the VINCI data reduction software (version 
3.0) by Kervella et al. (2004) employing the results based on 
the wavelets power spectral density. The calibration of the visi- 
bility values was performed as in Paper II using a weighted av- 
erage of the transfer function values obtained during the night. 



3.3. VLTI/VINCI results 

Table [3] shows the observational details together with the re- 
sulting calibrated squared visibility amplitudes for each series 
of y Sagittae interferograms. The listed errors include the scat- 
ter of the coherence factors of the single scans, the errors of the 
adopted diameter values of the calibration stars, and the vari- 
ation of the obtained transfer function during each observing 
night. 

As a first characterisation of the /T-band stellar angular di- 
ameter, we compute the equivalent UD (7 = yu°) and FDD 
(/ = fi l ) diameters, as for our NPOI data. The broad-band 
squared visibility amplitudes for the VINCI bandpass (X-band) 
are computed as 



IVatI- 



S~FlSl\V{v)?d V 
5^ n Sidy ' 



(4) 



where V(v) is the monochromatic visibility, F v is the stellar 
flux (assumed as Planck radiation) and S v the VINCI sensitiv- 
ity function including the transmission of the atmosphere, the 
optical fibers, the VINCI /^-band filter, and the detector quan- 
tum efficiency. 

Note that in the case of VLTI/VINCI the squared visibility 
amplitudes are integrated (Eq. |4}, while in the case of NPOI 
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Table 3. Details of our VLTI/VINCI observations of y Sagittae 
(date and time of observation, spatial frequency, azimuth angle 
of the projected baseline (E of N)), together with the measured 
squared visibility amplitudes and their errors. The last column 
denotes the number of successfully processed interferograms 
for each series. The effective wavelength for our y Sagittae ob- 
servations is ~2.19^/m. For each date of observation, we list 
the equivalent uniform disc (UD) diameter obtained from only 
the data of the specific night. Using all data together, we obtain 
an equivalent UD diameter of @ud = 5.93 + 0.02 mas, or an 
equivalent FDD diameter of ©fdd = 6.69 ± 0.02 mas. 



UT 


Sp. freq 
[1/"] 


az 
[deg] 


V 1 




# 




28 June 2002, UD 


= 5.91 ±0.03 mas 




05:16:57 


131.14 


136.63 


1.825e-01 


8.074e-03 


161 


05:22:40 


130.08 


136.65 


1.828e-01 


8.852e-03 


152 


05:32:34 


128.15 


136.73 


2.077e-01 


7.578e-03 


172 


06:44:22 


111.69 


139.42 


3.112e-01 


1.079e-02 


80 




8 July 2002, UD = 


= 5.89 ±0.04 mas 




05:12:10 


124.21 


137.05 


2.241e-01 


7.199e-03 


397 


05:20:13 


122.44 


137.27 


2.426e-01 


7.846e-03 


418 




11 July 2002, UD 


= 5.98 ±0.04 mas 




03:06:45 


142.16 


138.28 


1.064e-01 


1.256e-02 


55 


05:29:49 


117.47 


138.08 


2.495e-01 


8.775e-03 


69 


05:34:22 


116.38 


138.31 


2.660e-01 


9.929e-03 


207 


05:40:41 


114.83 


138.64 


2.878e-01 


8.524e-03 


297 


05:47:01 


113.25 


139.02 


2.977e-01 


1.228e-02 


198 




15 July 2002, UD 


= 5.94 ± 0.05 mas 




04:22:40 


128.76 


136.70 


1.926e-01 


6.824e-03 


90 


04:41:27 


124.90 


136.98 


2.223e-01 


1.315e-02 


139 




8 August 2002, UD 


= 5.92 ± 0.03 mas 




03:05:49 


125.17 


136.95 


2.179e-01 


7.678e-03 


469 


03:12:20 


123.77 


137.10 


2.244e-01 


9.005e-03 


436 


03:18:34 


122.39 


137.27 


2.309e-01 


7.979e-03 


468 


03:24:34 


121.03 


137.47 


2.471e-01 


1.901e-02 


178 


04:13:56 


108.94 


140.21 


3.376e-01 


1.212e-02 


423 




12 September 2002,0 


ud = 5.99 ± 


0. 1 1 mas 




01:26:22 


116.44 


138.29 


2.607e-01 


1.547e-02 


190 


01:47:17 


111.24 


139.54 


3.407e-01 


4.631e-02 


105 




18 September 2002, 


ud = 5.88± 


0.05 mas 




00:32:53 


123.38 


137.15 


2.347e-01 


1.081e-02 


428 


00:40:53 


121.59 


137.39 


2.393e-01 


1.177e-02 


407 


00:47:14 


120.14 


137.61 


2.667e-01 


1.751e-02 


210 



Table 4. Fit results of our VINCI data to UD and FDD models. 



Model 


Diameter 




xi 


UD 


UD = 5.93 ± 


0.02 mas 


0.63 


FDD 


FDD = 6.69 


± 0.02 mas 


0.63 



the visibility amplitudes have first to be integrated and squared 
thereafter (Eq.[3}, see the note in Sect. 12. 31 

Table @] lists the obtained diameter values for our VINCI 
data. Figure [8] shows our obtained VINCI squared visibility 
amplitudes of y Sge together with the best-fitting models with 
parameters listed below in Table|5] Since our VINCI data cover 
only one bandpass and only data of the first lobe of the visibility 
function, it is -contrary to our NPOI data- not feasible to con- 
strain the limb-darkening effect solely based on these VINCI 
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Fig. 8. Measured y Sagittae squared visibility amplitudes ob- 
tained with VLTI/VINCI in June to September 2002, together 
with the synthetic visibility curves of the best fitting models. 

data. This is also reflected by equal xl values obtained for UD 
and FDD models as well as by the virtually identical model 
visibility curves in Fig. [S] 

The increased equivalent UD diameter with respect to the 
shorter NPOI wavelengths is consistent with the general trend 
of decreasing strength of the limb-darkening effect with in- 
creasing wavelength. A detailed comparison of our data to 
model atmospheres follows below in Sect.|4] 

Analysis of calibration uncertainties In order to test and ver- 
ify the calibration uncertainties that are used in our analysis, 
we investigate the night-to-night variation of the obtained di- 
ameter values. All derived single nights' diameter values and 
uncertainties are consistent within 1 .3<x with the value obtained 
from all data together (5.93 mas) as well as with the weighted 
mean of the single nights' values (5.92 mas). 

This confirms that our diameter value is reliable and that 
our estimate of uncertainties is realistic. The obtained high- 
precision (0.3%) UD and FDD diameter values of ©ud = 
5.93 + 0.02 mas and Ofdd = 6.69 + 0.02 mas can thus be used 
without further uncertainties. 

Additional possible systematic errors that are constant over 
time scales larger than covered by our analysis, i.e. about 2 
months, can not be ruled out. Such systematic errors could in 
principle be related to the calibration of the interferometric ar- 
ray and the instrument, such as the calibration of the baseline 
length or the effective wavelength. Such uncertainties are not 
expected to represent a considerable source of error. 

4. Comparison to predictions by model 
atmospheres 

4.1. Employed model atmospheres 

We compare our measured visibility data to predictions by the- 
oretical model atmospheres in order to calibrate and test these 
models, and to derive fundamental stellar parameters of y Sge. 
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We use plane-parallel ATLAS 9 (Kurucz ll 9931 1 as well as 
plane -parallel and spherical PHOENIX (Hauschildt et al.[l999l 
model atmospheres to calculate synthetic visibility data, as 
done in Papers I&II. We refer to the descriptions in Papers 
I&II for more details on the employed model atmosphere files 
and their use. Differences between ATLAS 9 and PHOENIX mod- 
els include different opacity tables, a different sampling of the 
model wavelengths, a different sampling of the angles (jj. val- 
ues), and convective overshooting that is taken into account for 
ATLAS 9, but not for PHOENIX. 

The most important stellar input parameters for the plane- 
parallel models are effective temperature T e s and surface grav- 
ity log g, and for the spherical PHOENIX models in addition the 
mass M. We use solar chemical abundance, as appropriate for 
local cool giants. The values of T e g, logg, and M are already 
well constrained for y Sagittae, as outlined in the Introduction 
(Sect.G), namely T eff -3768 K, logg ~ 1.06, M ~ 1.3M . 

The closest model of the ATLAS 9 grid is the one for 
T e « =3750K and log g =1.0 (see Papers I&II for details on 
the model file used). We have constructed a corresponding 
plane-parallel PHOENIX model atmosphere with parameters 
T e s = 3750 K, \ogg =1, as well as a spherical PHOENIX model 
atmosphere with parameters T e ff =3750K, log g =1.0, and 
M = 13 (see Paper II for details on the model files). 

4.2. Calculation of synthetic visibility data 

We take into full account the bandpasses of our observations 
by integrating the synthetic visibility data of monochromatic 
intensity profiles for each spectral channel of NPOI and for 
the Zf-bandpass of VLTI/VINCI. This ensures that the synthetic 
visibility values fully resemble the true bandpasses used for the 
observations and that they include the model-predicted effects 
from atomic lines and molecular bands for each of our spectral 
channels. Monochromatic visibility values at frequency v are 
calculated as (cf. Davis et al. 2000 Eq. 6 from Paper I, Eq. 1 
from Paper II) 

V LD (v)= f (I v (ji)/I v (0))Jo(n®LD(B/X))fidfi. (5) 
Jo 

Here, I v (fi)/I v (0) is the normalised tabulated intensity profile, 
which is an output of the model atmosphere. Jo is the Bessel 
function of first kind and order 0; ©ld is the limb-darkened an- 
gular diameter at which the intensity profile reaches 0; B is the 
projected baseline length. Note that the evaluation of this inte- 
gral is vulnerable to numerical artifacts. We chose to use a lin- 
ear interpolation of the irregularly tabulated model values 
onto a regular grid of 1000 fi values between and 1 . The eval- 
uation of the integral was then performed using the Romberg 
method. Numerical results were checked against analytical re- 
sults for UD and FDD cases, and the resulting visibility func- 
tion for other cases was inspected for irregularities. 

Broad-band visibility values integrated over the bandpasses 
of our NPOI spectral channels and VINCI sensitivity function 
are calculated using Eqs.|3]and|4] respectively. 



4.3. Calculation of best fitting angular diameters 

We calculate the best fitting limb-darkened (0% intensity) an- 
gular diameter Old as described above for each of these three 
model atmospheres (plane-parallel ATLAS 9, plane-parallel 
PHOENIX, spherical PHOENIX models) and for each of our two 
data sets (NPOI and VLTI/VINCI). For each model fit, we treat 
Old as the only free parameter, and use all NPOI visibility data 
(squared visibility amplitudes, triple amplitudes, and closure 
phases), a total of 644 data points, simultaneously. The fit is a 
standard least-square fit, and optimises the total value of all 
644 NPOI data points. 

As discussed in Sect. 3.4 of Paper II, models based on 
plane-parallel geometry are optically thick from all viewing 
angles, the intensity steeply dropping to directly at the stel- 
lar limb. A plane-parallel model has, by definition, an atmo- 
sphere with an negligible thickness relative to the stellar ra- 
dius. Therefore, since any depth in such an atmosphere has a 
radius effectively equal to the stellar radius, a Rosseland diam- 
eter ©Ross in such a geometry is equivalent to the limb-darkened 
(or 0% intensity) diameter ©ld- 

Intensity profiles based on atmosphere models with spher- 
ical geometry, exhibit an inflection point and steepest decrease 
at radii smaller than the outermost model radius. The Rosseland 
mean optical depth increases slowly for increasing angles \i. 
Here, the ratio of the Rosseland diameter 0r oss and the 0% in- 
tensity diameter ©ld differs from unity, and this ratio Cr oss /ld 
is model-dependent and can be derived from the structure of the 
model atmosphere. This value depends on the definition of the 
outermost radius Rq of the model. Rq of the spherical PHOENIX 
model used here is given by the standard boundary conditions, 
which are a continuum optical depth of le-6 at 1.2 /mi and an 
outer gas pressure of le-4 dynes/cm 2 (see Paper II). 

4.4. Results and discussion 

Results Table |3 shows, separately for our NPOI and VINCI 
data sets, the resulting best-fitting angular diameter values 
based on the different considered model atmospheres, together 
with the corresponding^ 2 values. For the spherical model at- 
mospheres, the 0% intensity diameter ©ld is transformed to the 
Rosseland diameter 0r oss as described above. 

The corresponding synthetic visibility data are compared to 
the measured data in Figs. [3] to [7] for our NPOI data set and in 
Fig.|S]for our VLTI/VINCI data set. 

Best-fitting angular diameters Based on the plane-parallel 
ATLAS 9 model and our NPOI data, we reproduce the limb- 
darkened diameter L d = 6.18 + 0.06 mas from Paper I, despite 
the greater usable wavelength range and higher precision of the 
current NPOI data. The error includes an adopted 1 % system- 
atic error due to the NPOI wavelength calibration (the formal 
error is 0.004 mas). The VLTI/VINCI diameter for this model 
atmosphere of ©ld = 6.08 + 0.02 mas is not well consistent 
with the NPOI diameter (« 2 <x difference). 

The plane-parallel PHOENIX model leaves the near-infrared 
VLTI/VINCI diameter quasi unchanged (6.09 mas compared 
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Table 5. Results for the fit of ATLAS 9 and PHOENIX model atmospheres to our interferometric VLTI/VINCI and NPOI data sets 
of y Sagittae. 



Model atmosphere NPOI (526 nm to 852 nm) VLTI/VINCI (2190 nm) 



ATLAS 9, plane-parallel, T cS = 3750 K, log g =1.0 


©LD 


= 6.18 ±0.06 mas 


©LD 


= 6.08 


± 0.02 mas 






xl = 2.2 






Xl = 0.6 


PHOENIX, plane-parallel, T tff = 3750 K, logg = 1.0 


©LD 


= 6.11 ±0.06 mas 


©LD 


= 6.09 


± 0.02 mas 






Xl = 2.3 






Xl = 0.6 


PHOENIX, spherical, T eS = 3750 K, logg = 1.0, M = 1.3 Me 


©LD 


= 6.30 ± 0.06 mas 


©LD 


= 6.34 


± 0.02 mas 






Xl = 2.4 






Xl = 0.6 




©Ross 


= 6.02 ± 0.06 mas 


©Ross 


= 6.06 


± 0.02 mas 






©Ross = 6.06 


±0.02 


mas 





to 6.08 mas) with respect to the plane-parallel ATLAS model, 
while it results in a smaller (by « 1 cr) visual NPOI diam- 
eter compared to the ATLAS model. A comparison of these 
models' temperature structures reveals that the ATLAS model 
exhibits a steeper temperature gradient relative to the plane- 
parallel PHOENIX near Rosseland optical depth unity. This 
steeper gradient leads to stronger limb darkening at NPOI 
wavelengths and consequently a larger best angular diame- 
ter. The shallower temperature gradient of the plane-parallel 
PHOENIX model leads to a better agreement between the NPOI 
and VLTI/VINCI diameters. 

Finally, the spherical PHOENIX model leads to a Rosseland 
angular diameter of Or oss = 6.02 + 0.06 mas for the NPOI 
data set and Ross = 6.06 + 0.02 mas for the VLTI/VINCI 
data set. The larger best-fit diameters for the plane-parallel 
PHOENIX model compared to the spherical PHOENIX model ap- 
pears to be due to model geometry. The agreement of NPOI and 
VLTI/VINCI data sets within their 1 cr error bars gives confi- 
dence in both, the atmosphere models and the accuracy of the 
results from NPOI and VLTI/VINCI. The weighted mean of 
the NPOI and VLTI/VINCI results is Ross = 6.06 + 0.02 mas. 

Shape of the visibility function The measured and model- 
predicted visibility functions are generally consistent. 
However, the obtained reduced xl values for the NPOI data 
between 2.2 and 2.4 are above unity, as would be expected 
for a perfect match. This indicates differences at the 2cr level 
between observed visibility data and the model predictions. 

These differences are most evident in (1) a lower second 
maximum of the measured visibility function with respect to 
the model prediction on the EW baseline (Fig. [3j, and (2) a 
flattened measured visibility function with respect to the model 
predictions at the blue end on the CW baseline (Fig. |4j. It is 
not yet clear if and by how far these deviations of measured 
and synthetic visibility functions indicate different details of 
the limb-darkening effect at visual spectral channels, or if they 
are caused by additional calibration uncertainties that are not 
included in the error bars. In particular the flattening of the 
measured visibility function at the bluest spectral channels on 
the CW baseline can most likely be explained by additional 
calibration uncertainties of our NPOI data, as the instrumental 



transfer function for these data exhibits a drop which may not 
be fully compensated. 

The obtained diameter values in Tableware not affected by 
possible calibration uncertainties since the best-fitting diame- 
ter for any given model atmosphere is mostly constrained by 
the position of the first minimum and the global shape of the 
visibility curve. 

At optical wavelengths including all our NPOI spectral 
channels, TiO absorption bands are very important for the mod- 
elling of atmospheres of cool giants. It has been shown that the 
use of different line list combinations of TiO and H2O leads 
to significantly different model structures and spectra, in par- 
ticular in the optical where TiO bands are important (Allard 
et al. I2000I I. A possible explanation for differences between 
our visibility data and the model predictions could thus also 
be mismatching opacity tables for the TiO bands and/or a mis- 
matching spatial structure of the layers where TiO molecules 
reside. 

In order to estimate the effect of a lower model strength of 
the limb-darkening effect on the obtained diameter value, we 
used a spherical PHOENIX model with T e ff = 3500 K instead of 
our favourite model with T e ff = 3750 K (other parameters un- 
changed). The height of the second maximum of the visibility 
function at a wavelength of 600 nm is reduced from a 0.0085 
(see the lower right panel of Fig. [3 to « 0.0070. The obtained 
best-fitting diameter value for the NPOI data set changes from 
©Ross = 6.02 ± 0.06 mas to Or oss = 6.03 + 0.06 mas, and 
is unchanged for the VLTI/VINCI data set. This shows that an 
imperfect modelling of the strength of the limb-darkening ef- 
fect within our uncertainties does not have a significant effect 
on our obtained diameter values for y Sge. 

Model atmosphere fluxes Figure |9] shows the measured flux 
of y Sge from Alekseeva et al. 119971 in the wavelength range 
from 0.4-l.Oyum, i.e. covering the NPOI range used in this pa- 
per. Also shown are the predictions by the model atmospheres 
with parameters listed in Table|5] The limb-darkened 0% diam- 
eter values ©ld derived from the fit to the interferometric data 
were used to scale the model SEDs. The spectral resolution 
of the model SEDs is convolved to the resolution of the data 
used, i.e. to 10 nm. The model-predicted flux curves based on 
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Fig. 9. Flux of ySge from Alekseeva et al. (19971 in the 
wavelength range of our NPOI observations, compared to the 
model atmosphere predictions with model parameters listed in 
Table [5] The limb-darkened 0% diameter values ©ld derived 
from the fit to the interferometric data were used to scale the 
model SEDs. The spectral resolution of the model SEDs is con- 
volved to the resolution of the data used, i.e. to 10 nm 



the three considered models are well consistent with the gen- 
eral shape of the measured flux, while the detailed description 
of the spectral bands and features differs between the different 
models and the measured values. These differences can most 
likely be explained by the treatment of TiO absorption lines 
which are important for the visual wavelength range and dif- 



ficult to model (Allard et al. 2000), as mentioned in the para- 
graph above. 

Deviations from circular symmetry Differences between data 
and models are also observed for the closure phases (Fig. 0. 
The observed smooth variation of the closure phases from 
to n instead of the expected instantaneous flip may indicate 
a small deviation from spherical symmetry as already men- 
tioned in Paper I. As our maximum spatial resolution reaches 
2.7 mas, and the stellar disc has a size of 0r oss = 6.06 mas, 
the stellar disc is well resolved with 2.2 resolution elements 
across the disc. As a result, our data is sensitive to small de- 
viations of circular symmetry. Such a deviation can in princi- 
ple be caused by surface features such as spots, an asymmet- 
ric shape of the photosphere or of more extended (molecular) 
layers, or a faint unknown companion. Our target y Sge shows 
a relatively high photospheric pressure scale height of H?q = 
RgaJzslg * 0.006R*. Relatively large-scale (*0.06R*) sur- 
face inhomogeneities caused by convection could thus be ex- 
pected (cf. Freytag et al. 1997 1. 



Spherical versus plane-parallel model geometry The syn- 
thetic visibility data based on the same PHOENIX models that 
solely differ by spherical versus plane-parallel model geometry 



Table 6. Revised fundamental parameters of the M0 giant 
y Sagittae based on the analysis of this paper. For the details 
of the calculation, see the text. 



Parameter 



Value 



Rosseland angular diameter 
Rosseland linear radius 
Bolometric flux 
Effective temperature 
Luminosity 
Mass 

Surface gravity 



Ross = 6.06 ±0.02 mas 

R Ross = 55 ± 4R e 

/ boI = (2.57 ± 0.13) x 10-" W/m 2 

r eff = 3805 + 55 K 

log L/L e = 2.75+0.08 

M= 1.4 + 0.4M o 

logg = 1.1+0.2 



are virtually identical (Figs. I3I7I I. Thus, the model geometry 
has no noticeable effect on the shape of the visibility and can 
not be constrained by our visibility data. However, the derived 
angular diameter values differ by 1 .5% for the visual NPOI data 
and by 0.5% for the near-infrared VLTI/VINCI data. The spher- 
ical geometry allows us to more precisely define the stellar 
Rosseland radius with respect to the outermost model layer, and 
is thus more reliable than the angular diameter obtained from 
the plane-parallel model. As already noticed in Paper II for the 
M4 giant ifrPhe, the 0% intensity (LD) diameter based on a 
plane-parallel model seems to somewhat overestimate the stel- 
lar diameter with respect to the Rosseland diameter based on 
a spherical model. Our present results indicate a wavelength- 
dependent amount of this overestimation. 

5. Summary and conclusions 

We have compared our visual y Sagittae NPOI visibility 
data for 19 spectral channels with central wavelengths be- 
tween 526 nm to 852 nm as well as our near-infrared /T-band 
VLTI/VINCI visibility data with effective wavelength 2.19/mi 
to a plane-parallel ATLAS 9, a plane-parallel PHOENIX, and a 
spherical PHOENIX model atmosphere. The stellar parameters 
effective temperature T e e, surface gravity log g, and M of the 
model atmospheres used were fixed a-priori based on previous 
information on this star. 

The spherical geometry of the PHOENIX model enables us 
to precisely define the Rosseland radius of the star with re- 
spect to the outermost model layer and thus the 0% intensity 
diameter. This model leads to consistent Rosseland angular di- 
ameters for our NPOI and VLTI/VINCI data sets. This agree- 
ment increases the confidence in the model atmosphere predic- 
tions from optical to near-infrared wavelengths as well as in 
the calibration and accuracy of both interferometric facilities. 
In addition, the consistent angular diameter derived from our 
VLTI/VINCI data on a night-by-night basis over a total range 
of about 2 months increases confidence in the given calibration 
uncertainties. 

The Rosseland angular diameter of y Sagittae of 0r oss = 
6.06 + 0.02 mas, based on the comparison of our NPOI and 
VLTI/VINCI data to the spherical PHOENIX model, corre- 
sponds to a Rosseland linear radius of /?r oss = 55 + 4R Q , de- 
rived with the Hipparcos parallax of n — 11 .90 + 0.71 mas. The 
error of the Rosseland linear radius is dominated by the uncer- 
tainty of the parallax, not by the precision of our interferometric 
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measurement. With the bolometric flux ft, i = (2.57 + 0.13) x 
10~ 9 W/m 2 (Sect. ^ an d the Rosseland angular diameter, the 
effective temperature is constrained to T e $ = 3805 + 55 K. 
Again, the major contribution to this error originates from the 
uncertainty in /b i and not from our interferometric measure- 
ment. The Rosseland linear radius and r e ff result in a lumi- 
nosity of log L/L Q = 2.75 ± 0.08. Placing y Sagittae on the 
Hertzsprung Russel diagram using these values, and compar- 
ing to stellar evolutionary tracks by Girardi et al. (200(5) as 
in PapeiTI (Fig. 1 of Paper II) we can estimate a mass of 
M = 1.4+0.4M o , and thus a surface gravity of log g = 1.1+0.2. 
Table [6] summarises our revised values of ySagittae's funda- 
mental parameters. 

The closure phases show a smooth transition from to tt 
rather than a sharp flip, which could be due to a small devi- 
ation from circular symmetry of the well resolved stellar disc 
due to surface features such as spots, an asymmetric extended 
molecular layer, or a faint companion. 
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